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Abstract 

Definition of mixing fronts and turbulent mixing zone. Development of the Richtmyer- 
Meshkov instability leads to mixing of two substances separated by contact boundary. 
Let the axis z is perpendicular to a plane of boundary and z = corresponds to an 
unperturbed plane. The turbulent mixing zone is limited by front z = h + (t) > from 
heavy fluid side and by front z — h-(t) < from light fluid. Mixture of two fluids is 
confined between these two fronts (layer h_ < z < h + ). 

Initial conditions. Before pass of a shock wave both liquids were still [were at rest: 
v(x,y,z,t < 0) = 0] and y) was perturbed contact boundary. We consider 

linear (small amplitude) perturbation: |V?7i n j| <C 1. In this case the shock wave passes 
perturbations quickly: t s <C Xim/\VMni\i where t s is shock wave passage (sonic) time, Aj n j 
is a typical wavelength (or space scale) of perturbations, and \r]xini\ is typical velocity 
after shock wave passage. 

Transition and asymptotic stages. Universal asymptotics follows the transition stage 
at late time t ^> \/\f)\i n %\. Functions h±(t) asymptoticaly transform to power law depen- 
dences. Mixing is accompanied by cascade of enlarging of dominant scale of a turbulent 
mixing zone. Under condition ft = exact scaling laws are: h?P cx t 2 ^ and h?P oc t 1//3 , 
where \i is density ratio, 2D and 3D mark dimension of space. These scalings are calcu- 
lated in our work. The scaling exponents 2/5 and 1/3 are determined by the mechanism 
of redistribution of z— component of momentum from small into the large scales. 

For n 7^ the exponents 2/5 and 1/3 define the bottom limit of a range of possible 
values for both functions h+(t) and h-(t). Our 2D direct numerical simulations at different 
values [A ~ 1 show that d\nh±/dt w 2/5. 

This work has been supported by RBRF (grants 02-02-17499, 03-01-00700) 
and scientific schools (NSh-2045.2003.2, NSh- 70.2003.1). 
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# 1. Introduction 



The Richtmyer-Meshkov instability starts to develop after crossing (passage) by a 
shock wave of a contact boundary. In case of small on amplitude initial perturbation 
evolution in time proceeds in three stages (we shall designate them I, II and III). At short 
acoustic stage I 1 formation of a near boundary velocity field completes. This field appears 
due to interaction of contact boundary with a shock wave. Near boundary field defines 
the further hydrodynamical 2 motion. 

Duration of an acoustic stage is about 



Where \ ini is typical transverse scale of initial perturbations, c™ m is minimal from sonic 
velocities after shocks or after expansion wave. In case of usual equation of state this 
is sonic velocity in heavy gas because pressure is continuous at contact boundary and 



The expression (1) includes space scale \ ini (not amplitude r)\ ini ) as acoustic interaction 
of the curved contact boundary with the curved shock waves continues while shock waves 
will not depart from boundary on the large distances > \ n i. Interaction is carried out 
by means of a series of oblique sound waves which move in the expanding layer between 
shocks and boundary. This interaction creates vorticity which defines motion at stages II 
and III 3 . 

At a stage II near boundary velocity field decays at removal from boundary on dis- 
tances which exceed Aj n j. The law of reduction in velocity in depth from boundary depends 
on transversal correlations of perturbations. This law qualitatively differs for periodic, soli- 
tary and random perturbations 4 . In turn decay of velocity in depth determines the law of 
expansion of a zone of turbulent mixing at a stage III. The given work is devoted to the 
analysis of distribution of velocity near boundary. 

At stages II and III entropy and vortical modes on the one hand and acoustic modes 

J In case of small initial perturbations "hierarchy "of characteristic times for the stages 
I, II and III is following: tj <C tjj <C tin. 

2 In case of small initial perturbations hydrodynamical motion of substance at stages II 
and III occur to velocities small in comparison with sonic velocity (incompressible motion) . 
At stage I at presence of a shock wave with Mach number Ma significantly larger than 1 
compressibility effects are essential. 

3 Baroclinic generation of vorticity at the boundary and generation of vorticity in the 
bulk by curved shocks (Schilling et al., 2004; Wouchuk, 2001a, 2001b). In the case of weak 
shocks we can neglect bulk vorticity. 

4 Exponential decay for periodic perturbations and power laws for random and solitary 
perturbations at the same effective depth <~ X ini of near boundary velocity field. 



t s ~ \ni/ C" 
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with another - are decoupled 5 . Therefore an influence of compressibility becomes non- 
significant. Duration of sonic stage (1) is small in comparison with main time scale 

t\ini ~ \ni/ W\ ini . (2) 

The time (2) is necessary for development of perturbations up to nonlinear or "saturat- 
ed "amplitudes. Displacement of boundary during short time (1) is small in comparison 
with X ini . 

An estimate based on the Richtmyer formula gives 

~ At + k ini r)^ ini W unpertC B, (3) 

where At + is Atwood number after (sign plus) shock, ki n i = 2n/\ in i, i]^ ini is typical 
amplitudes of boundary perturbations (after shock) when we shift at distance ~ Aj„j in 
transverse direction, W unper tcB is unperturbed velocity of contact after shock, W unper tcB ~ 
c s . From equation (3) follows, that a smallness of t s (1) in comparison with t\ ini (2) it is 
caused by linearity of initial perturbations. 
Let 

Vxim ~ \ni (4) 

(nonlinear initial perturbations) , here r]^ ini is amplitude of surface perturbation before a 
shock. In this case we cannot separate stages I and II: t s ~ t\i n i, w\i n i ~ W unpert cB ~ 
c s . Hence, the stage III "will start "with high initial velocities making an appreciable 
share from sonic velocity. Nevertheless, received below asymptotical expansion laws of a 
turbulent mixing zone do not vary in case of large perturbation (4) in comparison with 
a case of linear perturbations. The matter is that eventually velocities are reduced with 
time, and become small in comparison with c s . 

Motion at a stage II is caused by movement of the primary vorticies ~ Aj„j created 
by a shock from surface perturbations at a sonic stage. As we will see, the essence of 
motion comes to light at introduction of concept about generations and terms primary, 
secondary, tertiary, etc. The question is the current generation of a chain of sign-variable 
vorticies of the size A curr (t) and about evolution of this generation. Such chain refers to 
also as dominating or large-scale structure. Evolution consists in growth of the size A curr (t) 
as a result of merging of vorticies. Each stage of merging connects parental generation 
and generation of first "descendants "of these "parents". Asymptotical law of expansion is 
established for far "descendants "of primary vorticies. Motion at a stage II has transitive 
character and is not universal. 

5 They are very different in velocities. Entropy and vortical modes are slow while acous- 
tic modes are fast. 
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Fundamental interest represents a finding of the law of expansion of a zone of mixing 
at large times t ^> t\ ini . This is an asymptotical turbulent stage (stage III) of devel- 
opment of the Richtmyer-Meshkov instability. Appears, present universal asymptotics, 
describing a stage III. Its representation makes the basic contents of the given work. As 
against laws of mixing at Rayleigh-Taylor and Kelvin-Helmholts turbulence when ex- 
ists simple self-similarity, here simple self-similarity is not present. In Rayleigh-Taylor 
and Kelvin-Helmholts cases there is a unique dimensional parameter: free fall accelera- 
tion g or a velocity difference Am accordingly. In the Richtmyer-Meshkov case there are 
two independent parameters Aj„j and u>Aim- Therefore simple self-similarity is absent. At 
such position presence universal asymptotics is rather not trivial bright phenomenon. The 
mechanism is below described, responsible for expansion of Richtmyer-Meshkov turbulent 
zone with self-similar hierarchy of generations. 

At stage III during each moment of time in a mixing layer the vortical structure 
dominates over the current scale A curr (t) ~ h(t), which grows with t (h is thickness 
of a layer of interpenetration). Evolution in time consists in change of one dominating 
structures by others by their enlarging (the inverse cascade, changing of generations). 
Externally it reminds a situation with Rayleigh-Taylor and Kelvin-Helmholts mixing. 
Actually mechanisms of enlarging is essentially different in Rayleigh-Taylor and Kelvin- 
Helmholts case on the one hand and Richtmyer-Meshkov with another. 

In the first case there is the plentiful external source which amplifyes subharmonics 
of the current dominating scale. In a Rayleigh-Taylor case it is gravitational energy, 
reserved in a heavy liquid. The heavy liquid is supported above light liquid due to pressure 
of light liquid. Therefore gravitational energy of a heavy liquid cannot be liberated in 
free falling. Downturn of the centre of gravity of system from heavy and light liquids 
occurs by exchanges or rearrangements of heavy and light liquids. Accordingly Rayleigh- 
Taylor instability is called also exchange or interchange instability. In Kelvin-Helmholts 
case a source is kinetic energy of motion which is external to a turbulent mixing zone. 
Alternation of generations consists in replacement of harmonics of the current dominating 
scale (parents) on subharmonics of these harmonics (descendants). 

In the second (Richtmyer-Meshkov) case an external source is not present. As it will be 
shown in the given work, fading Richtmyer-Meshkov motion is supported due to statistical 
redistribution of a longitudinal momentum from small scales in large scales. Thus in the 
second case (as against the first) energy of movement in the large scales is less than 
energy of movement, Which has been concentrated in small scales on the previous stages 
of development of turbulent mixing. 

In # 2 the statement of a problem is formulated. Connection between disintegration 
of break and Richtmyer-Meshkov instability is considered (generation of vorticies on con- 
tact at disintegration of break). It is spoken about "processing"of surface perturbation 
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r) into near boundary velocity non-uniform field v. The example of typical stochastic 
function with the finite scale of heterogeneity A in j is presented. Functions from this class 
characterize general initial perturbations. 

In # 3 characteristics of a mixing zone are described: (i) Average thickness of a mixing 
zone in which an interpenetration of phases takes place, (ii) large-scale structure (sequence 
of vortical dipoles), and (iii) small-scale details inside large-scale dipoles. 
In # 4 value of the pressure created in a liquid by the slowed down bubble is empha- 
sized. The important examples with a periodic chain of bubbles and a solitary bubble are 
considered. 

In # 5 the pressure and velocity fields near a random chain of bubbles is analysed. 

Comparison of periodic, solitary and random arrangements of bubbles is instructive. 

In # 6 spectral characteristics of a velocity field are calculated. It is found, that the major 

long wavelength part of a spectrum depends on wave number as a power law. 

In # 7 the velocity field arising directly after passage of a shock wave is analysed. 

In # 8 the law of expansion of a turbulent mixing zone into a heavy liquid is obtained 

(h + oc t 2/5 ). 

In # 9 dependence on time of large-scale velocity fluctuations is discussed (changeability 
with approximately constant step on hit in process of change of dominating structures in 
the inverse cascade). 

In # 10 the remark on mixing is made at arbitrary value of Atwood number. 
In # 11 stimulated and spontaneous turbulence are compared. In a stimulated case expan- 
sion of a turbulent mixing zone is connected to action of beforehand created large-scale 
perturbations. In a spontaneous case mixing develops from perturbations which have finite 
initial scale ~ \ ini . As against stimulated case it is a problem about decay of Richtmyer- 
Meshkov turbulence. It is interesting to consider motion on large times When thickness 
of a mixing layer is large in comparison with initial scale, h + \ ini . 
In # 12 received results are generalized on a three-dimensional case. The law 

h+ oc t 1/3 

of expansion of mixing front is found. 

# 2. Statement of problem 

In a problem about Richtmyer-Meshkov instability examine the movement arising 
owing to passage of a shock through contact boundary between two different gases. In 
fig. 1 straight lines 1 and 2 correspond to contact boundary (a straight line 1) and a shock 
wave (a straight line 2). 
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At the moment t — the wave reaches boundary. At t > from boundary in the 
different sides depart or two shock waves (a case 1), or a shock wave and a rarefaction 
(expansion) wave (a case 2) (see, for example, reviews in monographies Inogamov, Demi- 
anov, and Son 1999, Inogamov 1999 and in articles Wouchuk 2001a, 2001b). What from 
these two cases is carried out depends on the density ratio, thermodynamic properties of 
gases and Mach number of a shock. Usually in case of 1 shock goes from light gas into 
heavy gas ("from light to heavy"), and in case 2 - on the contrary. Thus the rarefaction 
wave leaves from contact boundary, moving on heavy gas. 

Problem about a shock wave, moving through plane contact boundary between two 
homogeneous semispaces, is a special case of the well-known important problem about 
disintegration of break. At the analysis of Richtmyer-Meshkov instability usually a case 
with substances, still before passing of a shock, is considered. Certainly, generalizations of 
this analysis are possible which include others cases of disintegration of break (generation 
of vorticity on and near the perturbed contact boundary at disintegration of general 
break). For example, it is interesting to consider a case in which before occurrence of 
shock two homogeneous semispaces move towards each other (Demchenko, 2004). 

Richtmyer-Meshkov instability arises when boundary between gases has perturbations. 
In Fig. 2 the example of the perturbed boundary is shown. The axis z (a longitudinal axis) 
is directed on a normal to unperturbed boundary. Below plane (2D) and three-dimensional 
(3D) spaces are considered. The coordinate x (x, y) is transverse coordinate in 2D (3D) 
space. 

Let's consider coordinate frame connected to unperturbed boundary. After passage 
of a shock through boundary there is a non-uniform near boundary velocity field v(f,t). 
It causes movings of boundary r](x, y, t) concerning position of unperturbed boundary 
(change of an amplitude of perturbations rj in time). In Fig. 3 movement of boundary in 
cases in a phase and an counterphase is shown (see, for example, Inogamov, Demianov, 
and Son 1999, Inogamov 1999). We compare here an angular phase of perturbations of 
boundary rj with a phase of velocity rj. Movement in a phase takes place at passage of a 
shock from light to heavy gas, and in an counterphase at passage of a shock from heavy 
to light gas (see Fig. 3 on which positive directions of longitudinal axes of coordinate z 
and velocity w are shown). 

In Fig. 4 velocity distribution v(x,z = 0,tfi x ) Or v(x,y = yfi X ,z = 0,tf ix ) along 
transverse coordinate x (or x, y) is shown in 2D (or 3D) case. As it was told, it is formed 
as a result of passage of a shock through the perturbed boundary. Plane z = corresponds 
to the current position of unperturbed boundary (dashed straight lines in Fig. 3). Velocity 
components u, v and w concern to axes x, y and z (compare Fig. 3 and 2). 

Because of the wavy form of the perturbed boundary (Fig. 2) with scale ~ Xmi the near 
boundary velocity field appears sign- variable with characteristic scales of wavelength \ ini 
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and amplitudes u^mi- Here the index "ini" specifies the initial moment of time. Let's con- 
sider a case of small perturbations vtini/^im ^ 1> when Richtmyer formula is applicable. 
The scale r\\ ini is a random mean square deviation of an interface at shift along boundary 
on distance ~ X ini . The index "plus "means that the given value is taken after passage of 
a shock. As it was spoken (# 1) found below power law expansion of a turbulent mixing 
zone remains fair for nonlinear initial perturbations (4). 

With dimensionless factor Richtmyer formula (3) looks like 

w n = FAt + k n r)+W unpertCB , w nm = FAt + k nm r)+ m W unpertCB . (5) 

It is used for recalculation of boundary perturbations rj into velocity field rj or w. In (5) F is 
dimensionless factor dependent in case of ideal gases from their adiabatic exponents, from 
temperatures of gases, and Mach number of a shock, see, for example, monographies Inog- 
amov, Demianov, and Son 1999, Inogamov 1999 with reviews of articles with calculation 
of factor F. In an incompressible case we have F — 1. This case corresponds to impulsive 
acceleration g oc 5(t). In (5) At + = , /i + is the density ratio after shock, W unpert cB is 
velocity of unperturbed boundary (a dashed straight line, Figs. 1-3) after shock in labora- 
tory frame in which substances before a shock are in rest. In case of small perturbations 
at a linear stage of development of instability a principle of linear superposition is valid. 
In (5) k n = 2nn/L, n is number of Fourier harmonic (2D), k nm = \J (k*) 2 + (km) 2 (3D), 
k nm is a wave number of a harmonic with x and y components k*, k y m in 3D geometry. 
We put our system into large (L ^> \ ini ) "box"with transversal sizes L (2D) Or L x L 
(3D). Perturbations continue periodically outside the end points of the "box". 

Comparison of amplitude of harmonic and a deviation due to shift. Values r/+ 
and r]+ m are Fourier amplitudes of harmonics. It is necessary to distinguish amplitudes r) n , 
w n of harmonics k n (k n = 2ir/X n ) from random-mean-square deviations due to transverse 
shift \ n 

V\n(t) = yj([ri(x,t) -rj(x + X n ,t)] 2 ) x , 

w Xn (z, t) = y/([w(x, z, t) - w(x + \ n , z, t)] 2 ) x , (f) x = J fdx/L. (6) 

Equations (6) can be generalized for 3D case. The random-mean-square deviation of ve- 
locity depends on velocity distribution which is formed by addition of many harmonics. 
Therefore for random functions / formed by an interference of harmonics f n with a ran- 
dom phase ip n we receive 

f\n ~ Vnfn- (7) 

# 3. Structure of turbulent mixing zone 
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Thickness of mixing layer. At t < there are two semispaces "h"and "l"in contact 
(Fig. 5a) filled with pure substances with density ph and p\. Thus a vertical (function of 
z) profile of average concentration c(z,t = 0) = (c) xy (see (6)) of light phase "l"is the 
"step"which is drawn by dashed straight lines in Fig. 5b. At moment of time t — the 
shock arrives at unperturbed contact boundary. After that in a picture of distribution of 
phases a new important zone appears - this is the turbulent mixing layer M (Fig. 5c), 
filled with mixture of substances "h"and "1". The quantity c(f,t) gives a volume fraction 
of light phase in an element of volume. Density of mixture which fills this volume is equal 

Pmix(r, t) = cpi + (1 - c)p h . 

After appearance of mixing zone the profile c(z,t) "is smeared" (Fig. 5b, a continuous 
curve) on thickness of a layer of mixing h = h + + h-. Coordinates z = h+ and z — h- 
(last points of interpenetration of phases) correspond to fronts of mixing: h + is a front of 
bubbles, and h- is front of jets (pi < p^). 

Transverse large-scale structures. Let us consider space structure of mixing zone. 
The one-dimensional profile of average concentration c(z,tf ix ) (Fig. 5b) reminds diffusive 
profiles. Intermolecular distances and lengths of free paths of molecules are small in com- 
parison with thickness h. Therefore molecular diffusive profiles are really one-dimensional. 
In a turbulent case magnitudes of the turbulent fluxes carrying phases along the z axis 
and rate of expansion of a mixing zone are defined by large-scale heterogeneities of mixing 
zone. It is important that the transverse sizes \ curr of these heterogeneities are large 

Acurr ~ h+. (8) 

Average flux (cw) xy of light phase "upwards" (to z > 0) and a heavy phase "downwards" ((1— 
c)w) xy are opposite. The mixing zone M (Fig. 5c) is two-stream region of counter motion of 
phases. Transport of "l"and "h"substances occurs by interchangement of volumes "l"and 
"h": h<->l (Fig. 6). Usually interchanges are connected with Rayleigh-Taylor instability 
where they free gravitational energy. Here (Richtmyer-Meshkov instability) interchanges 
(Fig. 6) occur because of inertia of motion, as the only possible form of motion. Certainly, 
any movings of liquid volume to an incompressible liquid are exchange motion. 

In space large-scale (Fig. 6) two counter stream motion is carried out through al- 
ternation of jets, "fingers"or columns of light "l"and heavy "h"liquids (Fig. 7a). This 
alternation causes changes of a sign on vertical velocity w as function of transverse or 
horizontal coordinates (Fig. 7b). Conformity of columns and signs on velocity is under- 
lined by the arrows connecting Figs. 7a and 7b. The correlator K cw (z, i) = ((c—l/2)w) xy ) 
is positive. The current thickness h of a mixing layer is defined on positions of last points 
h± (Fig. 7a). Horizontal arrows show designing these points on an axis z. 
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Sign alternation of velocity and alternation of columns are caused by a chain of vortical 
"spots" (Fig. 8a) in which vorticity turns off. Signs on circulation of "spots "in a chain 
alternate. It is shown in Fig. 8a with the help of the bent arrows in circles and in Fig. 8b 
by signs plus and a minus in circles. Vortical pairs (dipoles) involve in a zone of mixing 
pure "l"and "h"liquids (arrows "l"and "h"in Fig. 8a). 

Base element of large-scale structure is the vortical pair or a vortical dipole (Fig. 
8) of the size (8). The pair the next columns is connected to each dipole. Occurrence of 
vorticities occurs because of a variation of a sign on the momentum transferred by a shock 
to a near boundary layer of substance. In turn the variation of a momentum is caused by 
wagging about average the perturbed contact boundary r)~(x,y). Vortical circulation is 
connected to the angular momentum of liquid. 

Vortical motion at a stage III (# 1) is described in approximation of an incompressible 
liquid (conservation of volume). Thus overflows upwards should be compensated (Fig. 6) 
by returning of equal volume of a liquid downwards: J L J L w(x, y, z = 0,t)dxdy = 0. At 
random perturbations transverse correlations between base elements (between large-scale 
dipoles) fade on scale \ curr (8). That base elements appear under the order of size of one 
size along all infinite chain is caused by transversal uniformity on the average of initial 
perturbations rf~ . Therefore compensation of fluxes of volume J J wdxdy occurs mainly 
at size (8). Necessity of compensation results to interchanges (Fig. 6) and to alternation 
of columns (fig. 7). 

Previous stages and fractal structure. The instant spatial picture is above sub- 
mitted. Its basis is the large-scale chain of the dipoles set in scales A (8) m w rms , where 

w rms (t) = [[ L [ L [w(x,y,z = Q,t)} 2 dxdy/L 2 f 2 . (9) 
Jo Jo 

Amplitudes of fluctuations (9) is of order of random mean square w\ (6) at dipole size. 

Let this instant picture concerns to the moment of time P. The index % = \a.{t/tjj) (2) 
runs consecutive stages of development of the inverse cascade (sequence of enlarging of 
scale A*). At each enlarging the index i changes on size Ai ~ 1. In Fig. 9 two large-scale 
chains of dipoles A* at stages i\ and i 2 (i 2 — i\ >> 1) are shown. At a stage f 2 in a dipole 
a 2 there are almost all dipoles a\ from a stage f 1 got in an interval x a2 < x < x a2 + \ l2 . 
Thus, the dipole a 2 is the claster from the specified dipoles a±. 

In Fig. 10a dipoles of two consecutive generations i\ and i 2 = i\ + Ai are shown. In 
Fig. 10b three generations i±, i 2 = i\ + Ai and 23 = 22 + Ai are shown. Large-scale vorti- 
cal rotation strongly deforms and stretches small-scale details connected to the previous 
generations (Fig. 11). The analysis of numerical calculations shows that these details of 
current dominating structure are more likely relicts of previous generations than result of 
Kolmogorov (3D) or Kraichnan fragmentations. In such hierarchical way internal fractal 
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structure of large-scale formations forms. This fractal is responsible for rate of mixing of 
phases at a molecular level, see Redondo and Garzon, 2004. 

After separation molecular mixture remains as an intermediate layer in experiments on 
separation with change of a sign on acceleration, Kucherenko et al., 2004a. The contents of 
light phase in unit of volume of a mixture miscellaneous (varies from up to 1). Studying 
an intermediate layer it is possible to establish what share of the given layer having this 
or that contents of light phase. 

To estimate a share of molecular mixed mixture in total amount of mixture in tur- 
bulent zone it is possible by use of chemical indicators (Linden, Redondo, Youngs, 1994; 
Kucherenko et al., 2004b). Recent work (Chertkov, 2003) is devoted to studies of process 
of turbulent fragmentation during development of Rayleigh-Taylor mixing. 

# 4. Pressures and momentum, At= 1 

Let's describe calculation of an exponent 6 in a case \i = 0, At = 1 (mixing with 
"vacuum"). At a stage III flow tends to universal (self-similar) regime (# 1) in which h+ oc 
t 9 . Calculation 9 is based on the law of momentum conservation. Transfer of momentum to 
liquid is carried out due to ram pressure pv 2 of bubbles "imprinting"into liquid. Therefore 
creation of momentum by pressure is interesting for us. 

Pressure distribution. Periodical case. Let's consider a periodic chain of bubbles. 
Periodicity means that all bubbles in a chain are identical. Bubbles are symmetric con- 
cerning the vertical average line (boo), see Fig. 12b. The important point b is an apex 
of a bubble. In this figure the symmetric half of period of flow is shown. To receive all 
flow it is necessary to complete symmetrically to the left half of bubble up to the whole 
bubble. The full chain is obtained by periodic continuation of the whole bubble to the 
left and to the right 6 . 

Velocities of lagrangian particles on a straight line (boo) are directed upwards (to 
z > 0), arrows B in Fig. 12a and b, "imprint"of bubbles. These velocities (in the straight 
boo) decrease in due course (deceleration of bubbles, inertial imprint of bubbles). Hence, 
on the straight line (boo) the vector Vp is directed upwards. Distribution of pressure 
p(x — Xb, z,tfi x ) on an axis z is shown in Fig. 12a, Xb is coordinate of the bubble apex 
b. In case of an incompressible liquid the reference mark of pressure can be any function 
of time. In region occupied by "vacuum" (the region below the boundary of bubble bj, 
Fig. 12b) pressure is homogeneous: p(f,t) = pcBif). Let's put pcb = (a free choice of 
"calibration "of pressure). Then function p(z) in Fig. 12a begins from zero p(z = z^) = 0. 

6 About periodic solutions in 2D and 3D spaces see, for example, reviews Inogamov, 
Demianov, and Son 1999; Inogamov 1999; and works Inogamov and Oparin, 1999; Inog- 
amov et al., 2001a. Here for simplicity the example in 2D a case is presented. 
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Function p(z) in Fig. 12a grows with growths (as it is told, dp/dz > 0). Asymptotically 
at the large removals from roughnesses of boundary bj (z ^> A) heterogeneity of pressure 
fades - pressure becomes homogeneous: 

V -»■ Poo(*) > Pcb = 0, 

above A is the period of a chain of bubbles (Fig. 12b). It is shown in Fig. 12a and c. 

In a liquid the gradient Vp is finite on boundary bj between liquid and "vacuum". In 
jets j (liquid deformations aside "vacuum") the gradient |Vp| is small. Vertical straight 
line (joo) is a middle line of a jet. In Fig. 12b the family of instant isobars 0, 1, 2 and 
3 is presented. Pressure accrues with growth of number of an isobar. Zero isobar is the 
boundary bj. 

Initial kinetic energy of motion of liquid gradually passes in kinet ic energy of jets. The 
momentum 

/ Poodt 
Jo 

passes into a momentum of the substance which have been "thrown out"in jets j. Let's 
result an example. Let the condition with plane boundary and one harmonic of velocity 
was initial: 

z = r](x, t — 0) = 0, v = V</?, (p(x, z, t — 0) = — (w ini /k)e" kz cos kx, z > 0. 

Then initial vertical (w = ip z ) momentum Pf^, dx / °° wdz is equal to zero. At t = 
velocities of apex points b (wb, Xb = 0) and j (wj,Xj = ir/k) are equal each other ±w ini . 
At t — > oo absolute value of velocity Wj is approximately twice larger than Wi n i, Inogamov 
et al., 2001a. Integral 

rL roc 

I(t) = lim | / (dx/L) / wdz\ 

L^oo Jo Jr){x,t) 

(where function r](x, t) gives the current position of contact boundary) accrues in due 
course t and tends to limiting value 

roo 

/ Poo/phdt. 
Jo 

So, first, in a periodic case we have Poo > pcb (the case /x = is considered), 

Poo- Pcb ~ PhwI- (10) 

Second, heterogeneity of pressure are connected to roughnesses of boundary. 

Let's note that at /i ^ "hole "or depression in pressure is connected with compact 
region of vortical mixing (pressure in vortical region is lower than Poo)- Let's consider 
mixing near a horizontal plane z = (3D) or a horizontal straight line z = (2D). When 
plane mixing zone is present (Fig. 13a) then in distribution of pressure "valley" (2D) 
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extended along an axis x is formed. Average () x (6) instant profile of pressure p(z) is 
shown in Fig. 13b. Differences of pressure Ap 1 ^ = p 1 ^ — p min ~ Ph w l and Ap 1 ^ ~ p\w 2 
correspond to two semispaces z — > ±00 filled still pure "h"and "l"liquids with densities 
Phi Pi-, Pmin is pressure in a minimum. At a stage t ~ tjj (2) we have wj, ~ Wj ~ 
Accordingly Ap^ > Ap 1 ^ as p < 1, p = p t / p h . 

Pressure distribution. Solitary case. Let's consider a case when into liquid semis- 
pace (z > if) the single or solitary bubble (Fig. 14a) penetrates, 77 is contact boundary 
between liquid and " vacuum "bjoo. In Fig. 14a the right half of a bubble is drawn. This 
bubble is symmetric at transformation of a horizontal axis x — * —x. Flow is symmetric 
concerning a middle vertical straight line (boo), b is the apex of the bubble. Penetration 
of a bubble is accompanied by "squeezing"of liquid jets j. On a vertical straight line (boo) 
velocity (arrows w in Fig. 14b) is positive while acceleration of Lagrangian particles (ar- 
rows "acc", "dec") changes a sign: particles which are near to the bubble are decelerated, 
and distant - are accelerated. 

The profile of pressure with a maximum is shown in Fig. 14b. We have put Pb = Pcb = 

Where p& is pressure in the apex b. On infinity pressure in a liquid tends to the pressure 
in "vacuum"pcs- The system of isobars 0, 1, 2, 3, 4 is represented in Fig. 14c. The zero 
isobar is the contact bjoo. Pressure grows with growth of number of an isobar. 

# 5. Pressure and velocity distributions in a random case 

Let's consider a random chain (Fig. 15). Useful example of a random chain (together 
with numerical calculation its correlation and spectral characteristics) is presented in 
work Inogamov et al., 2001b. Horizontally homogeneous chains (which are made from 
comparable bubbles) are the most interesting. Let Aj n j is the average bubble size. The 
individual sizes of bubbles A, (Fig. 15) are randomly scattered about average value Aj„j 
with a dispersion AA ~ \ ini . Let's speak "periodicity about randomness": Q = \ ini /A\ ~ 

1 where Q is good quality (Q-factor). 

Penetration of horizontally unbound bubble chains raises pressure upon infinity p^ 
in comparison with pressure upon boundary pcb- Pressure p^ is created by ensemble of 
bubbles which is operating as the whole. This pressure is a measure of an average ram 
pressure pv 2 of bubbles in a chain. Concerning squeezing of infinitely removed liquid the 
action of the random chain is similar with periodic chain (# 4). 

In case of identical bubbles (a periodic chain) pressure in liquid is less than p^ (Fig. 
12, if z < 00) and heterogeneity of pressure is horizontally periodical. The pressure created 
by solitary bubble is more than p^ (Fig. 14). Movement of random ensemble creates a 
pressure field p(r,tfi X ) with a limit Poo{tfix) and random fluctuations in which pressure 
can be both or higher, or lower than limiting value Poo. 
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Important that fluctuations p form the hierarchical structure penetrating deeply into 
bulk of liquid in space. It is connected to formation of statistical groups of bubbles with 
different size of a ram pressure, average for the given group. 

"A groups"and "A fluctuations "of pressure. Let's describe hierarchical structure 
of heterogeneity of an instant field of pressure. Let's suppose that the ensemble is placed 
in "box"of the size L ^> \ ini with periodic boundary conditions at the box edges. Let's 
break length L into equal pieces A ^> X ini containing many bubbles ("A-group" of bubbles). 
Groups of bubbles on pieces A are various (Fig. 15). Influence of bubble group from a piece 
A reaches a depth of the order A from boundary into liquid. 

In Fig. 16a the figure 1 marks average position of contact boundary ("vacuum"is 
below). The length of the pieces 3 is proportional to pressure p\ in the center of a square 
A x A Worth the bottom side on the boundary. Changeability of this pressure is caused 
by dissimilarity of "A-groups"of bubbles. The braces 2 show average pressure {p) x on 
distance A/2 from the boundary. 

Creation of momentum of liquid in square cells A x A is connected to fluctuations of 
pressure 



as the homogeneous component (p) x does not make movement. Fluctuations of pressure 
(11) exponentially fade at distances from boundary exceeding A. Therefore at an estima- 
tion of momentum it is necessary to take cells with height of the order A. 
Let's write down an estimation of momentum 



In (12) A 2 is the area of a square, and w\ is an increment of velocity of substance in 
this square due to pressure action. At the left side of equation (12) there is a change 
of momentum in time per unit length in a transverse direction. On the right side in 
parentheses the average ram pressure of bubble is written down. Product of this pressure 
on \ ini gives average force per unit of transverse length. With such force bubble presses 
on liquid. Statistical multiplier as a root from an average number of bubbles in a "A- 
group"iVA gives estimate of random fluctuations of pressure (11), N\ = A/A^. The last 
multiplier on the right side of equation (12) is time of existence of bubble generation of 
scale Xini- 

As we see, momentum bubble groups of scale A written down in the right part of 
balance (12) diverges at large A because of the statistical factor oc y/\. It is caused 
by statistical accumulation of momentum with growth of the size of "A-groups". Thus 
velocity of motion w\ oc I\/ S\ decreases, here I\ and S\ are momentum and the area of 
a "A-group". 



Ap A = Px - (p) 



(11) 




(12) 
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The sign of velocity fluctuations W\ (12) randomly varies from one square A x A to 
another (Fig. 16b). Signs w\ and Ap\ (11) coincide. 

Hierarchy of " A-groups"and dependence from z. Above we have considered one 
near boundary layer of thickness ~ A in which the random velocity field w\ is created 
by " A-groups"from ~ N\ bubbles. Horizontal (random changeability) and vertical (expo- 
nential fade) scales of spatial heterogeneity of this field is of order of the length of pieces 
A, Figs. 15 and 16b. 

To find velocity field from a random chain (Fig. 15) it is necessary to consider set 
of splittings into pieces A of different scale. The hierarchy of splittings (or hierarchy of 
"A-groups") looks like: Ao, Ai, A2..., Where 

A n = fl n Aj„j, (13) 

a > 1, a ~ 1, for example, a = 2. With each scale (set by the index n) particular near 
boundary layer is connected. On the one hand scales with different n are independent 
(there are not correlations between them), on the other hand the final velocity field w(z) 
(from a chain) develops of imposing all layers n. 

For an example we shall consider "A-groups "from ten (on the average) bubbles (N\ = 
A/Aj ni = 10). Changeability on coordinate x values A^a (11) and w\ oc Aj>a (12) is 
generated by non-equivalence of neighboring tens, Fig. 16a and b. This changeability 
fades on depth ~ 10Aj ni . According to the formula (12), the amplitude of the given 
changeability makes w\q oc 1/10 3 / 2 . 

Let's compare "A-groups "from tens and hundreds bubbles. The amplitude of change- 
ability of velocity at shifts on coordinate x on 100A in j is less (w W q oc 1/100 3 / 2 ). Therefore 
on depth ~ 10Aj n j fluctuations W100 are imperceptible on a background of fluctuations W\q 
(closer over contact is dominated with fineer splittings). But on distances from boundary 
~ 100Aj„i fluctuations w\q exponentially decrease, and to dominate begin fluctuations 

WlOO- 

Told it is illustrated in Fig. 17 on which the spatial map of vertical velocity w on a 
plane of movement (x, z) is shown. Three consecutive scales n — 1, n, n + 1 hierarchies (13) 
(small, average and big squares) are represented. Velocities w n -i,w n and w n+ \ dominate 
over different distances. To negative velocities there correspond the shaded squares. Thus, 
the field w(x, z) looks like randomly branching tree with a trunk L, sharing on a (13) 
branches which in turn share for a branches everyone, etc. 

Dependence of amplitudes of fluctuations of velocity on depth w(z) is formed as result 
of "interference "of splittings A n . Total action of hierarchy of splittings (13) can be written 
down as 

W ( Z ) =J2 W n ex P K J2 A « 3/2 eX P K ' W ^ ~ Wxini ~~fk " 

n V \i ' n V J Z 1 Z 
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Step-type behaviour or lacunarity (13) (lacuna or gap) for hierarchy was necessary 
for simplification of explanations. It is possible to consider an interference of splittings 
with fine step of a geometrical progression a — 1 + e, e <C 1. Velocity field is smooth; 
its particular shape is defined by arrangement of bubbles in a chain (Fig. 15), instead 
of hierarchy (13). Thus of "stain "or "spots "of negative velocity w have smooth contours 
(not squares, as in Fig. 17). In Fig. 17 the diagram of splitting of consecutive scales of 
these "spots", instead of their form is submitted. 

The Fig. 18 supplements the formula (14). Points n — 1, n and n+1 mark the centers of 
squares, compare Figs. 17 and 18a. These centers are on distances A n /2 from the boundary 
submitted by a horizontal piece. The specified points are transferred on a plane (w, z) 
(Fig. 18b) with the help of three horizontal arrows connecting the left and right parts 
Fig. 18. Coordinates of these points are (w n oc l/\f/ 2 ,z n = A n /2). Through sequence of 
the transferred points the curve w(z) (14) is drawn. As we see, small-scale fluctuations 
w\ (12) have the big amplitude whereas large-scale will penetrate deeply into thickness 
of liquid. 

Now it is possible to compare instant velocity fields w(x,z) from a periodic chain (# 
4), from a solitary bubble (# 4), and from a random chain (# 5). These qualitatively 
various three cases definitely supplement each other. In the first case decay on z occurs 
exponentially. Power law dependence from z is connected to a solitary (or lonely) bubble, 
however distribution of velocity on coordinate x is non-uniform. Velocity field from a 
random chain, first, is homogeneous on transverse distances large in comparison with 
the size of a bubble, and, second, as a power law depends on coordinate z. As it will be 
visible below, power law dependence w from z results in the power law of expansion of a 
turbulent mixing zone (h oc t e ). In a periodic case mixing quickly fades (displacement of 
bubbles h oc hat). 

# 6. Spectrum 

Let's calculate Fourier spectrum of an instant velocity field w(x,z = 0) created by a 
random chain. The scale 

Xcurr? kcurr 27f/ \curr? ^Icurr Xcurr 

divides a spectrum on long-wave (n < n curr ) and short-wave (n > n curr ) parts. The size 
\ curr corresponds to scale of dominating structure of the current generation of large-scale 
vorticies or the dipoles developed by the moment of time t ~ t curr , see Fig. 7, 8 and # 
3. For questions examined in the given work the short-wave part is unimportant. The 
long-wave part is determined by fluctuations w\ (12) at A > A curr . 

The scale of change of velocity w\ at length A, see (6), is connected to amplitudes of 
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harmonics w n by equation (7). Let's write down w\ ~ \^nw n , where numbers of harmonics 
n correspond to scale A : n ~ ri\ = Lj 'A. From here we receive 

w n = w Xini (\ ini /L) 3/2 £ n cos i/j n n, (15) 

where £ n ~ 1 is dimensionless number and ip n is the random phase. It turns out, that a 
long-wave spectrum is linear under number n. At t ~ t curr the formula (15) is valid up to 
n < n curr 7 . At a derivation was considered, that velocity field w(x,z,t ~ ti n i) keeps the 
kind (under the order of size) in scales A > A CMrr down to the moments of time t ~ t curr 8 

# 7. Random chain and a shock 

We research velocity field w + (x, z) formed by a shock during short sound time t s 
(1) near a random chain of linear perturbations rf (# 1, # 2). Equation (5) connects 
harmonics of function r] + (x) and harmonics of initial velocity 9 Sequence of three steps 
of calculation following. First, we define the Fourier spectrum 77+ of random surface 
perturbations rj + (x) w . Second, under Richtmyer formula the amplitude of harmonics w n 
is defined. Thirdly, on a spectrum w n is function w + (x, z) defined and dependences of 
this function from x and z are analyzed. Corresponding calculations are executed in work 
Inogamov, 2002. It is shown, first, similarly (15) spectrum w n is linear on n, second, 
dependence w from z is given by the formula (14). 

# 8. Decay of velocity with depth and expansion law 

Let's address to a question on evolution of velocity field w(x, z, t) (the given section 
and the following). After a push connected to a shock wave, there is velocity field w + (x, z) 

7). We admit, in the further field w + remains "frozen"in liquid. In the following section 
we shall return to this assumption. Let's explain, that means "frozen". Basis is that (as 
it was spoken # 5) small-scale modes have large velocities, and large-scale will penetrate 
deeply into liquid. 

The field of velocity w + (x, z) with enlarging cells (Figs. 17, 18a) and the power law 
decay on z (Fig. 18b) is formed. 

At such position the assumption about "frozen"cells is natural. Really, turning time 
of large formations is large 

t x ~ \/w x oc t 5/2 . (16) 

7 Instead of for n < riini = L/Xi n i. In time t ~ t curr the maximum of a spectrum is 
displaced from n ~ rii n i in n <~ n curr . 

8 To discussion of this question are devoted # 8 and # 9 below. 
9 Initial for a stage of movement in a mode of an incompressible liquid. 
10 It appears, the long-wave part of a spectrum rj n does not depend from n, r\ n oc n°. 
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Change of velocity w\ on depth z ~ A occurs after penetration of contact boundary on 
this depth h + ~ A. If it so it is possible to write down 



^~ w (,)U + ~i/4 /2 . 



From here follows the required law 11 of decay of Richtmyer-Meshkov turbulence 



h + ~wV 5 ni \^t e , e = 2/5. (17) 

Large-scale formations enter game one by one in the hierarchical order (13) and with the 
big delay (16). This gradual inclusion of more and more large cells defines the law h+(t) 
(17). 

# 9. Accumulation of momentum and logarithmic corrections 

Velocity field w + (x, z) is formed (# 7) for a time interval ~ t s (1). Further in time 
~ t\ini (2) "balls" (Fig. 8) primary vorticies are turned off - the dominating structure of 
the first generation is formed. At t ~ t s displacement of contact boundary are small (still 
there is no turning boundary). The random chain of sign- variable vorticies Aj n j, Wxim forms 
velocity field w(x, z,t ~ t\i n i) in near boundary layer. Important that fields w + (x, z) and 
w(x, Zjtxini) under the order of size are identical. 

Further process of change of one generations by the following develops. Let's con- 
sider the current generation with scale of dominating structure \ cur r- The random chain 
Kurr, wxcurr forms velocity field 

w(x,Z,t ~ txcurr) ~ W X (t Xcurr ) ~ W Acurr (A CMrr /A) 3/2 (18) 

in surpassing scales A > A curr . Substituting in (18) following from the law (17) estimations 

, /, % dh+(t\ curr ) 

■^curr ~ " + \t\curr ) 5 U^Xcurr ~ ^ j 

11 Several phenomenological approaches was earlier offered (K-model and buoyancy-drag 
model) for the approximate estimation of an exponent 6. Usually reduction of this exponent 
in comparison with value 2/3 (conservation of energy) is connected with viscous losses of 
energy on small scales. The review of these approaches is given in work Inogamov 2002. 

In work of Schwarz et al., 2000 the exponent 9 for 2D and 3D cases has been extimated 
with the help of bubble envelop model. The envelop model essentially distinguishes from 
K-model and buoyancy-drag model. Values 82D ~ 0.4 and 9zd ~ 0.2 have been received 
by the envelop model (Schwarz et al., 2000). 

As against the listed approaches in the present work the exponents 02D (17) and 6*313 (# 
12 below) are calculated precisely. The mechanism of expansion of mixing zone is opened. 

Appeared, dynamics of expansion is defined by long wavelength (A > A CMrr ) fluctuations 
of velocity spectrum. At the moment of time t curr these fluctuations are at a linear stage 
of development (displacement of liquid particles 77a on scale A are small in comparison 
with A). 
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we express velocities w(x, z,t ~ t\ cu „) or W\(t\ curr ) in (18) through parameters of primary 
vorticies \ n i, W\ ini . It turns out, that in the rest by the moment t ~ t\ curr undeveloped (or 
"unsaturated") 12 scales A > A CMrr all generations of dominating structures create under 
the order of size the same field of velocity ~ w + \ 

At discussion of velocities w(x, z, t\ curT ) or w\(t\ curr ) the question was amplitudes of 
fluctuations. Let's consider now the phase information. On amplitude velocities created 
by different generations are of the same order of magnitude value. It is asked, whether the 
fields created by different generations are identical on a phase? In other words, whether 
the shaded cells on a map of velocities from different generations (Fig. 17) coincide? 

There are three opportunities. First, the large-scale part (A > A CMrr ) fields w(x, z) can 
practically not vary in time. In this case the law of mixing (17) is carried out. Second, 
structure of shaded cells formed by different generations can coincide under the order of 
magnitude. Thirdly, this structure can randomly vary from generation to generation. 

Let's consider last two opportunities. Let's recollect, that w\(t\ curr ) is an increment 
of velocity connected to an increment of momentum for a time interval ~ t curr . Means 
for the finding of velocity field developing by the moment of time t ~ t curr it is necessary 
to summarize all increments caused by activity of the previous generations. 

In case the second opportunity is realized, we receive 

W\(t Xc urr) ~ hl(A curr /% ni ) W + . (19) 

In the third case we have 

W\{t X curr) ~ \j ^{Karr / \ni) W + . (20) 

Accumulation of increments results in logarithmic amplification. The logarithmic correc- 
tions are connected with a little bit strengthening of mixing. In view of it instead of the 
law (17) it is received 

/i + oc(lnt) 2 / 5 t 2 / 5 , h+ oc (hit) 1 / 5 t 5 / 2 

in cases (19) and (20) accordingly. 
In 3D case we have 

/i+ocQnt) 1 / 3 * 1 / 3 , /i+ocOnt) 1 / 6 ^ 3 . 

# 10. Influence of inertia of light phase 

In the previous sections the case fi = (At = 1, contact to "vacuum", [i is density 
ratio) was investigated. The case \i ^ (the account of density of light phase) is difficult, 
12 The latent or "frozen"scales. 
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and demands separate consideration. Position at p ^ in 2D geometry is difficult. In this 
case kinetic energy of an incompressible liquid conserves 13 for the case of small viscosity 
(Re ^> 1). In a case p — the viscous heat production can be neglected at large values 
of Reynolds number Re (conservation of kinetic energy). However difficulty with excess 
of energy in small scales 14 at p = does not arise. The matter is that the fly of fast 
jets in "vacuum"carries away surplus of energy. At such position promotion of front h + 
aside heavy liquid is defined by the law of conservation of momentum which results in 
the formula (17). 

Flow with the power law of expansion h + oc t e is self-similar. Thus in self-similar 
variables 

t°n, t'-^wn 

spectrum (# 6) is steady-state 15 . The degree scaling amplitudes, appears from an equation: 
y/nw n ~ h + at (7) and # 6 (definition of n curr ). 

In self-similar flow the profile of average density p = (p)j_ is function not of two (z, t) 
but of one variable 16 , p(z/h + ) or p(z/t e ); the mark _L means averaging on transverse 
coordinates. The coordinate z is counted from initial position of boundary (# 1, # 2). 
At p = function p quickly decreases aside "vacuum". It becomes very small already at 
value of dimensionless coordinate (— z)/h + (t) about several units. 

First portions of the substance (which have been thrown out by primary generation 
of bubbles) flying ahead aside "vacuum"move under the law h- ~ —w\i n it. Distribution 
of density near to front /i_ is not described by the self-similar solution p(z/h + ). 

#11. Stimulated and spontaneous turbulence 

Addition of initial long wavelength spectrum. Problem I is above considered in 
which the long wavelength part (A > A curr ) of a spectrum w n or w\ is formed owing to 
an interference of random small-scale modes. Let's address to other problem (a problem 
II) 17 . In it to a random initial chain of surface perturbations \ n i, ri\ ini the long wavelength 
spectrum is added 

^/3+3/2 

Vn W = 1/3+1/2 C0S ^" n/3 ' ^<n<n ini = L/\ ini , (21) 

13 Let's note, that in a case [i = viscous dissipation at Re— > oo is absent in 3D geometry. 

14 Formally the law of conservation of energy results in value 9 = 2/3 which essentially 
exceeds the value (17). 

1B It is steady-state on the average. Thus there are fluctuations of amplitudes about 
average values and fluctuations of phases. 

16 Certainly, self-similar asymptotics is established in process of growth of number of the 
generations replaced after generation of primary vorticies Xi n i , w\i n i . 

17 Below problems I and II refer to as problems about stimulated and spontaneous tur- 
bulence. 
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where (£„) = 1, ip n is a random phase, L is length of "box"(# 2). Let's write down also 

V LW = ( ^ + 3/2 /jL /m/2 )n/3+ l/ 2) e LW = V LW /X „ ( A ^ L )/?+3/2 ^+3/2 _ (^/^/W (22) 

where is a relative amplitude of surface perturbation. It characterizes an angle of a 
deviation of a tangent to perturbed surface from unperturbed position. 

At A ~ A* we have r)% w ~ A (22) - the assumption of linearity of perturbations is 
broken. As we are interested with long wavelengths A > A*, should be 

(3 > -3/2. (23) 

Value 

(3 = -3/2 (24) 
corresponds to a homogeneous spectrum of perturbations 

/ n - . e hom L/n 3 / 2 oc l/\k ± \ 3 / 2 , ^ w oc 1/n, e{ w ~ e hom = const, (25) 



in 



where k± is transverse or horizontal wave number. Homogeneous noise (25) are dimen- 
sionless in the sense that do not contain parameter A* with dimension of length as in this 
case it is carried out r)% w ~ A. They are homogeneous in space of wave numbers as have 
no the allocated scale (from here the name). These noise play the important role in the 
theory of the turbulence connected to Rayleigh-Taylor instability. 

Let's "sew"a long wavelength part of the initial spectrum set by formulas (21), (22) 
with perturbations "located "at wave number axis X ini , w\ ini 18 . Then amplitudes w^ w 
(curves 1 and 3 in Fig. 19) and u^i™ (a curve 4) are of the same order at A ~ Aj„j. 
Let's write out a condition of sewing together e{ w ~ r\\mil\ni having taken e^ w in this 
condition at A ~ Aj„j. Let's express parameter A* through Aj„j with the help of the given 
condition 

A* ~ (^Ami/Aj„i) 1/(/3+3/2) Xini- 

We see, at r\\mi < Xini for the spectra limited to a condition (23) is carried out A* < Aj„j, 
see Fig. 19. 

Excluding A* it is received 

x /3+3/2 x /3+3/2 / > \ 13+3/2 

n LW _ \ni B LW _ \ni /H-l/2 LW _ ^ini ] „/3+3/2 ( ^ 

'In ~ t mi ^/3+l/2 ' 'A ~ ml £8+1/2 ' A ~~ t%m 1 £ j n ' l zu / 



18 The short wavelength tail A < Xini is not significant for dynamics of system at t » 
txini- The long wavelength tail created by "located"at k ~ kini perturbation (Fig. 15) 
looks like w n (x n (15). 

Long wavelength spectra (15) and (21) are different. 
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Where the notation ej n j = rj\i n i/\ ini is entered. Equations (26) describe added 19 long 
wavelength initial "noise". They replace expressions (21) and (22). 

Development of Richtmyer-Meshkov turbulence at presence of added long 
wavelength "noise". In Fig. 19 perturbation \ n i,w\ini ( a curve 4) "located"on an axis 
n is shown. Large-scale fluctuations of the "located "perturbation generate (## 5-9) a 
long wavelength part of a spectrum. It is represented by a curve 2. Thus, the spectrum 
representing association of curves 2 and 4 is connected to the "located"perturbation. 

In Fig. 19 added long wavelength "noise "is represented with curves 1 and 3 concerning 
to different values of an index (5 (26). For an example from a curve 1 we have (5 < 
0, and for a curve 3 it is carried out (5 > 0. At (3 = dependences e\ from A are 
identical to added long wavelength "noise" (a curve 2) and for large-scale fluctuations of 
the "located "perturbation (a curve 2). Let's note, that under condition of (23) for these 
dependences the limiting equation e\ — > is carried out at A — > oo. 

Interest represents a case — 3/2 < j3 < when added long wavelength "noise "surpasses 
large-scale fluctuations of the "located"perturbation. Such "strong"noise defines dynamics 
of mixing at t > t\ ini . 

The harmonics of velocity corresponding to this "noise "are defined by Richtmyer 
formula (5) 

w™ ~ (2n/L)nri w W unpertCB oc n 1+/3 . 
For large-scale fluctuations of velocity it is received 

w x ~ y/nw n oc r/ +3/2 . (27) 

Let's note, that according to Richtmyer formula it is carried out W\ oc e\. 

Let's write down dh + /dt ~ w\ curr (see # 8 and # 9). We have h + ~ A CMrr (Figs. 6-8). 
From here we receive 

h ~ 1/(W2) x (/3+3/2)/(/3+5/2) l/(/3+5/2) .l/(/3+5/2) ^ 

~ A CMrr ~ € ini A ini W unpertCB 1 ~ 

~ A£ +3/2)/(/3+5/2) ^l{i? +5/2) t 1/(/3+5/2) , (28) 

as wxcurr oc 1/Af^ /2 according to (27). 

Homogeneous noise (24), (25) is the extreme case of the large-scale noise stimulating 
mixing. In this case during the initial moment we have W\ ~ w\i n i at anyone A > Aj„j. 
Hence, the formula (28) becomes h + ~ w\ ini t. 

Comparison of stimulated and spontaneous mixing. There is a certain analogy 
between classes of turbulent movements at Rayleigh-Taylor and Richtmyer-Meshkov in- 
stabilities. There is a division on stimulated and spontaneous classes in Rayleigh-Taylor 

19 Added in the initial data to "located"on wave number k initial perturbation. This is 
the combined ("located"and long wavelength) or compound perturbation. 
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case (Inogamov, Demianov, and Son, 1999; Inogamov, 1999) and analogue of this division 
in Richtmyer-Meshkov case 20 . Let's remind sense of concepts about stimulated and spon- 
taneous motions. At spontaneous turbulence unstable boundary is initially "seeded"with 
the perturbations "located "on k or A near X ini . Are interested in an self-similar regime 
which is established asymptotically at h + \ ini . Flow tends to a self-similar expansion 
regime of a mixing zone 

h± = a ± Mgt 2 (29) 

with factors a± not dependent on a concrete shape of "located "on k initial perturbations 
("forgeting"of initial perturbations ~ Xini)- 

We call stimulated a regime with long wavelength stimulation. Imposing at t = long 
wavelength power law spectra such as (26) introduces scale A* and thus "breaks "self- 
similar dependence (29). The special place belongs to homogeneous "noise" with pertur- 
bations of boundary 77 a kind (25) and/or velocity perturbations of a kind 

™n ~ € v hom ^L/n oc l/\k ± \, w x ~el om jA^\. (30) 

In 3D a case (see also #12 below) we have the following expressions for homogeneous 
noise 



thomL . ._ 2 thom-L v \J AtgL 1-3/2 / c 9 nV 

Vnm ~ -T72— OC \k±\ , T] X ~ — , W nm ~ € hom /2 OC 1 . (8.30) 

^nm -< v nm Nnrn 

Homogeneous noise do not contain dimensional parameter A*. Hence, self-similarity (29) 
it is kept. However, with "noise" (25), (30), (30)' dimensionless numbers or amplitudes 
are connected e;, om , e^. In such situation factors a± become functions a±(ehom, e v hom ) of 
these dimensionless amplitudes describing intensity of "noise "stimulating mixing. 

Real perturbation is a combination of spontaneous and stimulated 21 perturbations. 
Dependence ot±(ehom, e nom ) is slowly growing function of the arguments (a curve "st"in 
Fig. 20). The factor af +st (e hom ) (a curve "sp+st") is limited from bottom 22 by the con- 
stant value a+ (a straight line "sp"). From the data on factor a+ suggested in report 
Youngs 2004 in Fig. 20 value of 0.03 is accepted. There is a value of dimensionless or 
relative amplitude a ^ which influence of stimulation becomes appreciable. 

20 According to it let's carry problem I (see the beginning # 11) to spontaneous 
Richtmyer-Meshkov turbulence and a problem II - to stimulated. 

21 In cases with small values of stabilizing factors (viscosity, heat conductivity, diffusion) 
spontaneous turbulence gets ineradicable character while stimulated turbulence - it is 
supervised by dimensionless amplitudes e^om- 

22 Due to reduction of thorn it is impossible to achieve that mixing will develope more 
slowly than spontaneous mixing (the spontaneous regime "props up "the combined modes 
"from below"). Similar situation exists for the Richtmyer-Meshkov mixing, see Fig. 21 
below. 
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It is very important that numerically the value e^ m is rather small. It was supposed 
in the previous estimations (Inogamov, Demianov, and Son, 1999; Inogamov, 1999) and 
has been confirmed in recent detailed numerical calculations (Youngs 2004). Smallness 
of the factor a+ and smallness of the threshold amplitude ^om strongly raise a role of 
a regime with stimulation. The amplitude e^ om gives value of deviations of a tangent to 
boundary on all scales as at homogeneous noise e\ does not depend from A. Smallness of 
^hom means that hidden or latent widezone perturbations are significant. It is caused by 
that these weak widezone noise can increase essentially the major characteristic of mixing 
- the factor a + . 

In reports Youngs 2004 and Dimonte et al. 2004 it is assumed, that available diver- 
gences between numerical (a + ~ 0.03) and experimental (a + ~ 0.05) estimations a + are 
caused by the latent stimulation of mixing in experiment. Such stimulation difficultly to 
supervise because of smallness of ^hom- The origin of eh om can be connected to vibration 
of experimental installations. The interesting program (Dimonte et al., 2004) for search 
of weak widezone noise in experiment is offered. 

Analogue of stimulation of mixing in Richtmyer-Meshkov case. At Richtmyer- 
Meshkov turbulence large-scale perturbations (26) intensify mixing. In stimulated regime 
(fig. 21) the exponent 9 is more, than in a spontaneous regime: 9 st > 9 sp , 9 st = (j3+5/2)~ 1 , 
see (28) 23 , 9 sp = 2/5, see (17), # 8, # 9. Borders of an interval of values of an index 
f3 form points j3 = and f3 = —3/2 on an axis j3. In Fig. 21 more to the right of the 
specified interval the forbidden area, see (23) settles down. More to the left - the area 
lays, in which stimulation it is possible to neglect (compare spectra 2 and 3 in Fig. 19). 
The point (3 = —3/2 (24) represents homogeneous or "dimensionlessnoise"(25). To them 
there corresponds the maximal value of function 9 st ((3.) 

The most important 24 is the point (3 — 0. It concerns to a mode which we shall 
conditionally name "spontaneous". The reserve here contains in the following. The com- 
parative analysis of Rayleigh-Taylor and Richtmyer-Meshkov problems is carried out. 
At Rayleigh-Taylor turbulence in a spontaneous case there is no dynamically essential 
long-wave tail existing on "all times"t > (to it is connected "forgeting"small-scale ini- 
tial perturbations). On the contrary, in the stimulated case such tail is present (constant 
presence of a part A > A curr of an initial spectrum). In this respect at Richtmyer-Meshkov 
mixing business is in another way. Namely, both in stimulated and in spontaneous cases 
there is constantly existing long wavelength tail determining dynamics. Hence, and in a 

23 In report Youngs 2004 stimulation of Richtmyer-Meshkov mixing is numerically in- 
vestigated. Power law initial surface perturbations of a kind (26) has been used. Change 
of the exponent is revealed at change of an index (3. 

24 As this point represents the most widespread variant of initial conditions With "lo- 
calization "of perturbations on wave number k (a random chain Xini)- 
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spontaneous case the certain "trace"of the initial data is kept. 

Richtmyer-Meshkov turbulence in a regime (17) is generated by random "located"on 
k perturbation (## 5-8). Such perturbation generates spontaneous mixing in Rayleigh- 
Taylor case. The basis is connected to it to name spontaneous a regime of mixing (17) in 
Richtmyer-Meshkov case. 



# 12. 3D case 



Let's generalize the description of turbulent mixing received above on 3D space. The 
horizontal plane with coordinates x, y (Fig. 2) is transverse plane. Let's limit a part of 
this plane by big (L ^> Aj n j) square L x L ("box", # 2). Let's consider a random two- 
dimensional lattice of perturbations on a plane x, y. Such lattice (3D) replaces a random 
chain (2D, # 5). Let the average transversal size of initial surface perturbations 

z = rj ini (x,y) (31) 

in this lattice is equal to \ ini and randomness is characterized by good quality Q ~ 1 (# 
5). 

Let's expand periodic in a square L x L function (31) in Fourier series 

oo oo 

Vini(%? U) = ijInmC nxC m y + Vnm^nxSmy + f] nm SnxC m y + Vnm^nxS m y) , (32) 
n=l m=l 

Where C nx = cos(2imx/L), S nx = sm(27inx/L). Amplitudes of Fourier harmonics in 
series (32) are calculated by integration 

ri nm = (4:/L 2 ) / dxdyrj ini (x : y)C nx C my . (33) 
jo Jo 

Similar formulas define other amplitudes which are included in a series (32). 

Long wavelength (n, m N ini = L/\ ini ) asymptotics of spectra r] nm looks like 

Vnm ~ inm cos ^ nm (4/ L 2 )r] rms \f ni ^J A 2 / \? ni yJ L 2 / A 2 ~ r] rms X ini / L ~ r) rms /N ini , (34) 

where (£ nm ) = 1, is a random phase, r\ 2 Tms = (1/ L 2 ) J L J L dx dyr]f ni , A = 2n/k, 
k 2 = (^«) 2 + (^m) 2 ; Ki = 27rn/L, k v m = 2nm/L are components of a wave vector on 
a plane (x,y). This vector corresponds to a harmonic (n, m). Asymptotics (34) follows 
from a statistical estimation of integral (33) of stochastic function (31). Occurrence of 
statistical factors is connected to stochastic fluctuations of this function \J\ 2 /\ 2 ni and 



L 2 /\ 2 . The estimation (34) is valid on an interval of values n, m from 1 up to ~ N, 
In 3D a case instead of (7) we have r]\ ~ \/n 2 + m 2 r] nm ~ i] rms \ ini / A, where n ~ m 
L/X. 



ini- 
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Let perturbation (31) is linear (small amplitude): i] rms <C X ini . Let's estimate velocity 
field after passage of a shock wave. We use Richtmyer formula connecting amplitudes 
of harmonics of surface perturbations r] nm and vertical velocity w nm . In 3D geometry we 
have 

w nm = FAt+(27r/L)v / n 2 + m 2 7]+ m W unpert cB, (8.35) 

compare with the formula (5). 

Let's find velocity of movement of primary bubbles wa™- Such bubbles have transverse 
scale ~ Xini. Let's substitute in the formula (35) expression for amplitude r\ nm (34) at 
n ~ m ~ N ini . In result it is received 

w nm (n, m ~ N ini ) ~ FAt + (r] rms /L)W U n P ertcB- (36) 

Let's write down 

Wrms ~ Wxini ~ W X {X ~ Xini) ~ N in iWnm{n, Ul ~ N ini ) (37) 

and 

w A (A ~ X nm > A ini ) ~ (L/X nm )w nm , Km = L / \/ n 2 + m 2 . (38) 

In the formula (38) size w\ is determined by the formula (6) generalized on 3D a case. 
Formulas (37), (38) replace the formula (7) in 3D a case. Fourier amplitudes w n (2D) or 
w nm (3D) are connected to the corresponding size w\ by means of the statistical factor of 
order of a root square from number of modes of scale A. In 2D case this number is L/X 
and in 3D a case - is (L/X) 2 . 

From formulas (36), (37) it is possible to define velocity of primary bubbles w\ ini . The 
estimation of this velocity looks like 

~ FAt + Nini(r] rms /L)W un pertCB- (39) 

Asymptotics of 3D Richtmyer- Meshkov mixing at late times is defined by a long 
wavelength part (n, m <^ N ini ) of a spectrum of velocity w nm . Let's find asymptotics of a 
spectrum w nm with the help of asymptotics of a spectrum r\ n m (34) , Richtmyer formula 
(35) and auxiliary expression (39). We have 



w nm ~ {Vn 2 + m 2 /N?Jw Xi ni- (40) 

It is useful to compare 3D expression (40) to the formula (15) valid in 2D case. 

From formulas (38) and (40) the asymptotics (A ^> Aj n j) of velocity fluctuations of 
scale A can be found. It may be written as 

^2 _|_ ^2 / A ' ' \ ^ 

W X (X ~ X nm > \ni) ~ "'Aim ~ ( ) W Xin i- (41) 



im 
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The hierarchy (13) of fluctuations (41) of different scales forms enlarging cell structure 
of velocity field w(x,y,z) similar to presented in Fig. 17. Three-dimensional cells have 
the sizes ~ A x A x A. As well as in 2D a case, decrease of amplitude of fluctuations of 
3D fields w(x,y,z) at removal from boundary is described by the power law. To tell the 
truth, an exponent is another. From the formula (41) it is received 

w(z) ~ Wxini{Kni/ zf , (42) 

compare with (14). 

Decay of velocity w(z) defines the law of penetration of mixing front h + (t). Let's write 
down dh + /dt ~ w(z ~ h + ) (see # 8). Substituting here the formula (42) and integrating 
we receive 

K ~ w^XfSt 9 , e = 1/3. (43) 

Let's compare laws (43) and (17) concerning to 2D and 3D to cases. We see, in a three- 
dimensional case Richtmyer-Meshkov turbulence fades faster and the front h + {t) moves 
ahead essentially more slowly 

where values 62D and O^d are given by formulas (17) and (43) accordingly. 
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Figures 
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Fig. 1. Unperturbed flow. It is shown approach of a shock wave (a straight line 2 and an 
arrow) to contact boundary (a straight line 1). 




Fig. 2. Richtmyer-Meshkov instability arises from surface perturbations (a continuous 
curve). Unpertubed position of boundary is shown by dashed straight line. Definition of 

longitudinal and transverse coordinates. 




Fig. 3. Baroclinic "processing"by a shock wave of surface perturbations r] into velocity 
field v. The form (continuous curves) and directions of movement (arrows) of contact 
boundary after passage of a shock for cases in a phase (at the left) and in a counterphase 
(on the right). Shaped straight line is position of unperturbed boundary. 
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Fig. 4. The most interesting example (by virtue of the typicalness) of the general 
stochastic function set in scales Aj„j and w\ini (w is a component of velocity on an axis 
z). The sign- variable field of velocity in the near boundary layer arising after passage of 
a shock through the perturbed boundary is presented. 




Fig. 5. The average distribution of light "l"and heavy "h"phases (a) prior to the 
beginning of mixing and (c) after the beginning of mixing. As a result of interpenetration 
of phases intermediate layer "M"with a mixture of "l"and "h"phases is formed. 
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Fig. 6. Involving into turbulent movement of volumes of light "l"and heavy "h"liquids. 




Fig. 7. (a) The instant picture alternating "l"and "h"columns. (b) Sign variability of 
velocity w(x,yfi X ,z = 0,tfi X ). The plane z = corresponds to position of unperturbed 

contact. 
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Fig. 8. Large-scale structure and a vortical chain, (a) "Blowing"of "l"and "h"columns 
by pairs of vortical spots, (b) A chain of large-scale sign- variable vorticities or a chain of 

dipoles A x h from pairs of vorticities. 




Fig. 9. Formation of a claster dipole from small dipoles of the previous stages. 
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Fig. 10. Compound character of large-scale dipoles at late 



stages of mixing. 






Fig. 12. Squeezing (p^ > 0) of bulk (+00) of liquids by ram pressure of bubbles pressed 
into it. (a) p(z) on a straight line (boo); (b) Increase of p on isobars (0 is a boundary 
of a bubble "bj"), 1, 2 and 3. At boundary "bj"pcB = 0. (c) p(z) on (joo). 





(a) 




Fig. 13. (a) The layer with the turbulent motion concentrated near a plane z = 0. (b) A 
profile p(z) with "hole"(/z 7^ 0) in the region of movement. 
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Fig. 14. (a) Penetration of bubble. Geometry of flow, (b) Distribution p lengthways 

(boo), (c) An instant field of pressure. 
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Fig. 15. Ensemble of bubbles A, with the average size Aj n j and groups of bubbles, got on 

pieces A (A-groups). 



1< 



3 * 



< 



u. 



^ — »»[*S — 



Co) 




(*) 

Fig. 16. (a) Changeability of pressure (pieces 3), created by "A-groups"of bubbles, (b) 
Random alternation of squares A x A with positive and negative (are shaded) increments 

of velocity w\. 
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Fig. 17. Hierarchy of changes of a sign of velocity w(x, z) : small-scale changeability 
near surface, large-scale - it is far in a bulk of liquid. 




Fig. 18. (a) Scales n — l,n, n + 1 and the squares connected to them, (b) "Designing"the 
centers of squares from the left figure on right. 
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Fig. 19. The combined perturbations from (1) added long wavelength "noise" (a curve 1) 
and (2) "located"on A a spectrum (a curve 4). Perturbations 4 generate a large-scale tail 

(a curve 2). 
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